Interplay between interferences and electron-electron interactions 

in epitaxial graphene 

B. Jouault^'^, B. Jabakhanji^, N. Camara^'^, W. Desrat^, C. Consejo^, J. Camassel^'^ 
Universite Montpellier 2, Groupe Etude des Semiconducteurs, 
cc074, pi- Eugene Bataillon, 34095 Montpellier cedex 5, France 
2 CNRS, UMR 5650, cc074, pi. Eugene Bataillon, 
34095 Montpellier cedex 5, France and 
^ CNM, Campus UAB, Bellaterra (Barcelona), Spain 

Abstract 

We separate localization and interaction effects in epitaxial graphene devices grown on the C-face 
of an 8-°off 4H-SiC substrate by analyzing the low temperature conductivities. Weak localization 
and antilocalization are extracted at low magnetic fields, after elimination of a geometric mag- 
netoresistance and subtraction of the magnetic field dependent Drude conductivity. The electron 
electron interaction correction is extracted at higher magnetic fields, where localization effects dis- 
appear. Both phenomena are weak but sizable and of the same order of magnitude. If compared to 
graphene on silicon dioxide, electron electron interaction on epitaxial graphene are not significantly 
reduced by the larger dielectric constant of the SiC substrate. 
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I. INTRODUCTION 



Graphene-based devices are exciting candidates for future generations of microelectronic 
devices. A promising technique to produce graphene at an industrial scale is the epitax- 
ial graphene growth from a SiC substrate, because these SiC substrates can be patterned 
using standard lithography methods. Thanks to recent technical improvements, the most 
delicate and intrinsic features of graphene, those reflecting the chiral nature of the quasi- 
particles, as, for instance, the so-called 'Half Integer quantum Hall effect-"- and the weak 
antilocalization^, have been recently reported for epitaxial graphene. 

In graphene, depending on the relative magnitude of intervalley scattering time and 
phase coherence time, either weak localization (WL) or weak antilocalization (WAL) has 
been predicted^. The phase interference correction to the resistance depends on the nature 
of the disorder. For epitaxial graphene samples, elastic scattering favorable for WAL can 
be caused by remote charges like ionized impurities in the substrate. On the other hand, 
atomically sharp disorder (local defects, edges) causes intervalley scattering and gives rise 
to WL. 

Experimentally, in two-dimensional gases, WL is often mixed with electron-electron in- 
teraction^ (EEI). As WL, EEI gives a correction to the Drude conductivity with a InT 
dependence. It follows that the experimental extraction and separation of EEI and WAL 
contributions are usually difficult^. For graphene, in the diffusive regime, EEI is expected 
to give this usual temperature dependence correction proportional to ln(T). EEI is also 
expected to be sensitive to the different kinds of disorders^. 

In this work, we take advantage of simultaneous measurements of the longitudinal and 
transverse resistances to invert the resistivity tensor. We can then separate the different 
mechanisms which give corrections to the Drude conductivity: i) a geometric contribution, 
which experimentally appears as a constant term in the longitudinal conductivity; ii) WL 
and WAL, by a comparison between the experimental magnetoconductance and the Drude 
magnetoconductivity, iii) EEI, whose temperature dependence can be analyzed in the mag- 
netic field range for which WAL has disappeared. 
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II. EXPERIMENTAL DETAILS 



We have used large and homogeneous single epitaxial graphene layers grown on the C- 
face of an insulating 8° off-axis 4H-SiC substrate. The graphene sheets have an elongated 
triangular shape of which quality and homogeneity can be easily checked using micro- Raman 
spectroscopy. On few selected samples, Cr/Au ohmic contacts were deposited to define 
Hall-bars with a rough geometry (see inset in Fig.l for details). Then magneto-transport 
measurements were done using a 14T magnet in a cryostat and a Variable Temperature 
Inset operated down to 1.5K. On the best samples, at high magnetic field, the half-integer 
quantum Hall effect could be observed up to the last plateau in the temperature range 1.5 
to 40K. For details, see Ref.— . In this work we focus on two moderately doped samples 
(5*1 and 52) with dimensions, carrier concentration, mobility and scattering times reported 
in table [H We work at low injection currents (lOnA- 1/iA), mainly at low magnetic fields 
{\B\ < 3T), low temperatures (1.5 K-200K) and we focus on the magnetic field dependence 
of the longitudinal and transverse elements of the resistivity tensor. 



III. BACKGROUND CONSIDERATIONS 



An overview of the results obtained for sample SI is shown in Fig.l. The experimental 
longitudinal and transverse magnetoresistances Rxx{B) and Rxy{B) for sample 5*1 are pre- 
sented in Fig. la and lb respectively, at different temperatures between 1.6K and 48K. Let 
us consider first. Fig. la. At low fields [B < 0.1 T), a negative magnetoresistance peak 
centered &i B = OT is observed in Rxx- This peak is typical of weak localization (WL). 
Indeed, this peak a linear dependence of Rxx vs ln(T) with a slope of the order of h/e^ (see 
inset of Fig. la) and the amplitude of the peak is ^pxx ~ 2p^^e^//i {pxx is the longitudinal 
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SI 
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5 1.1 5000 0.57 0.1 


0.065 


0.03 


52 


40 


10 0.8 11000 0.48 0.15 


0.11 


0.057 



TABLE I: for sample 51 and 52: length L between lateral probes, width W (in /-tm), hole con- 
centration (in lO^^cm"^), mobility (in cm^V~"'^s"^), scattering times T^r, Tir and Ttr (in ps) and 
diffusion constant (in m^s""^). 
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FIG. 1: (color online) (a) Rxx{B) at different temperatures for sample SI. All curves crosses at 
B ~ 2T and have a behavior. Inset shows the temperature dependence of Rxx at B= T 
(black open squares) and B= 0.45 T (solid red squares), (b) Rxy{B) at different temperatures. 
The slope of the Hall voltage increases when T decreases due to electron-electron interactions. 
Inset: photograph of sample SI, edges of the graphene layer have been indicated by a dashed line. 

resistivity). At higher magnetic fields, reproducible fluctuations in the conductance G are 
also observed, with amplitudes AG ~ 0.5e^//i for 5*1 and AG ~ O.le'^/h for 5*2. 

At magnetic fields \B\ < 0.4T, a smooth depression is observed in Rxx, which we at- 
tribute to weak antilocalization (WAL). Because this WAL is barely visible in Fig. la, 
we also present additional experimental results in Appendix A. Experimentally, WAL is 
blurred because fluctuations of conductance are also present, and also because the positive 
magnetoresistance of the WAL is superposed to another positive magnetoresistance with a 
pronounced parabolic dependence, which is well observed on the whole magnetic field range 
\B\ < 3T presented in Fig. la. The theory of weak localization is based on the diffusion 
approximation and does not hold for magnetic fields much higher than Btr= h/AeDrtr ~ 
100 mT, where Ttr is the transport relaxation time and D the diffusion coefficient. This 
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means that the parabohc background above Btr should not be attributed to WL features. 
It must have a different origin. Since the device is being far from having the shape of an 
ideal Hall-bar (see inset in Fig. lb) we ascribe this parabolic magnetoresistance component 
to magnetic deflection of the current lines. This is because the lateral probes are invasive 
and because the graphene layer under the lateral contacts is likely to have different mobility 
and carrier concentration^i. For details, see Annex B. 

Beyond geometric and interference corrections, a third correction to the resistivity mani- 
fests. It comes from electron-electron interactions and shows as: i) the persistence of a ln(T) 
linear temperature dependence of Rxx, in a magnetic field range where interference effects 
are suppressed (see inset of Fig. la for a field of -B = 0.45T); ii) a variation of the Hall slope 
(Fig. lb) as a function of T, which cannot be explained by variations of the carrier density, 
as the hole gas is strongly degenerate [Ep/kBT > 20 for T < 50K, where Ep is the Fermi 
energy); iii) a crossing of all Rxx{B) resistances taken at different temperatures aX B ^ 
(Fig. la)i^. 

The justification of the last two points, originally given in^, is as follows. EEI give no 
correction to the transverse conductance (Jxy. It gives a small correction 5a^^ to the longitu- 
dinal conductance a^x- This correction does not depend on B if the magnetic field is smaller 
than a critical field Bs < TiksT /2^b- This condition if fulfilled for all the temperature and 
magnetic field ranges of this study. Because of this correction, EEI gives rise to a negative 
magnetoresistance in the first order in da"^^ ^ axx'- 

P.. ~ 1/ao - (1 - ^^^B'')5a''/al (1) 

and a variation of the Hall slope ^Pxyjpxy ~ —Ibo^^ja^ where ctq is the conductivity at 
i? = 0. These relations are derived by inverting the conductivity tensor: 

'^xx Icy. 

^ XX ' xy 

^ XX ' xy 

It has been established recently^ that the separation of EEI and WL^ is simplified if 
the magnetoconductivities axy{B) and (Jxx{B) are used, rather than the resistivities. As 
already stressed, the conductivity of the graphene layer is different below the Cr/Au contacts 
and, strictly speaking, inverting the resistivity tensor is incorrect because the device is 
inhomogeneous. However, experimentally, the geometric correction is small, appears only 



in the longitudinal magnetoresistance as a parabolic correction, and consequently geometric 
corrections to the longitudinal conductivity appear as a -B-independent shift 5a^ . 

The magnetoconductivities for sample SI are plotted in Fig. [2] at different temperatures. 
The conductivities have been obtained from Eq. [2] and |3l where the resistivities have been 
estimated by p,, ^ {W/ L){R,,{B) + R,,{-B))/2 and p,y ^ (R^yiB) - R,y{-B))/2. The 
lower inset of Fig.2 shows that all a^y^B) taken at different temperatures collapse on the same 
curve, as expected for EEL This also confirms that both mobility and carrier concentration 
are constant over the whole temperature range and, as the geometric correction only depends 
on T via the Hall angle fiB, that Scr'^ does not depend on the temperature. Finally, the T- 
independence of the transverse conductivity indicates an effective separation of EEI {6a'^'^) 
and interference {6a^^) corrections over almost the whole 5-range. At low fields \B\ < 
Btr, a T-dependence of a^y due to WL is expected but is not observed, being beyond our 
experimental resolution. 

The large decrease of a^xiB) between and 2T in Fig. [2] is due to the magnetic field 
dependence of the Drude conductivity = ao/(l+/i^-B^). We express the total conductivity 
by 

= a^ + da"" + 6a^ + Sa""^. (4) 

When WL is negligible: \B\ > 0.5 ^ Btr, we impose 6a'^'' = and Eq. |U for a given 
temperature, simplifies + C, where C is a constant which incorporates both 

geometric and EEI corrections. This last formula gives very good fits to the conductivity 
at all temperatures. As an example, the blue dotted line in Fig. 2 is the fit for the T=5.6K 
data in the interval 0.5-2 T, where ctq, /i and C are the fit parameters. Transverse and 
longitudinal resistivities can be fitted separately but still give very similar mobilities and 
concentrations, which are reported in Table I. 

The upper inset in Fig. 2 is an enlargement of the low-i? data at T = b.QK and evidences 
the weak antilocalization peak which takes place around the weak localization centered at 
5 = 0. A similar weak antilocalization can be detected for all temperatures. The WL 
theory does not take into account the modification of the current paths due to magnetic 
field. In other words, before examining the WL contribution to the conductivity, we should 
first subtract the Drude fits + C (given by the dotted line in Fig. |2]) to the observed 
longitudinal conductivities^^. Experimental results for samples 5*2 and 5*1 are shown in 
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Fig. [3l for different temperatures. 



IV. CORRECTION DUE TO WEAK LOCALIZATION AND ANTILOCALIZA- 
TION 



In the diffusive regime in graphene, 5a^^ is given by^: 

Txh 



+ 2Bsr J V + ^sr + Bir 

•1 I 1 



(5) 



where F{z) = In z + + -), ^1/ is the digamma function, B^^sr,ir = JT'/4:DeT^^sr,ir and T^^sr,ir 
are the coherence time, the intervalley scattering time and the intravalley scattering time 
respectively. For simphcity, we identify intervalley to short range (sr) scattering and in- 
travalley to long-range (Ir) scattering. We take constant scattering times Tsr and Tir over the 
whole temperature range and only is allowed to depend on T. Neglecting the warping^^, 
we also impose r^-^ = -|- r^"^. We then fit the conductivities at different temperatures by 
Eq. |5l Results are indicated by dashed lines in Fig. [3] (for sake of clarity, fits at high tem- 
peratures have not been reported). Best fit is obtained with the values reported in Table I. 



We estimate the short-range scattering length L''' = ^/Dt^ ^ 100-140 nm. The long-range 



scattering length L'^ = yDr^ 70-90 nm is even shorter. These lengths are comparable to 
what is found in exfoliated graphene on Si02 substrate^^. They are also comparable to the 
distance between the SiC steps below the graphene layer. It is often stated that exfoliated 
graphene is much more disordered at the edge of the SiC steps^^. Therefore step edges could 
be the main source of scatterings in these samples. 

The phase coherence time obtained from the fit is plotted in Fig. HI Apart from a 
saturation at low T, is roughly proportional to 1/T between 10K-50K with a slope equal 
to ^ 50 ps K for samples 5*1 and 5*2. We conclude that t^^ obeys the usual temperature 
dependence for electron electron scattering in the diffusive regime: 

T-^=PkBT\ng/hg, (6) 

where g is the reduced conductivity: g = (j^h/e^ and the empirical coefficient /3 is 1 for 5*1 
and 1.4 for 5*2. Similar observations have already been done both for epitaxial^ and exfoliated 
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FIG. 2: (color online) longitudinal conductivity Gxx versus B at different temperatures for sample 
SI. Blue dashed line is the best fit of the conductivity at T= 5.6K, on a i?-interval 0.5-2 T. The 
fit is done according to cr^xi^) + ^ (^^^ text). The upper inset is the enlargement of the area 
around B ~ OT. The lower inset shows that all transverse conductivities (Txy{B), taken at different 
temperatures, collapse on a single curve. 
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FIG. 3: (color online) (a) Corrections to the conductivities da^'' + C at different temperatures for 
sample SI. The T-dependent vertical shift C has been kept for clarity. WL fits according to Eq. [5] 
are also indicated, (b) Similar analysis for sample 52. 

graphene^*^. The coherence length = >^/Dt^ is readily calculated and plotted in fig. 4b. 
At low temperature, slightly exceeds l/xm. 



V. CORRECTION DUE TO ELECTRON ELECTRON INTERACTION 

We attribute the temperature dependence of the conductivity to EEL The Drude mobility 
is constant over the whole temperature range, therefore there is no temperature dependence 
of the conductivity induced by electron-phonon scattering. This is not surprising, as the 
temperature dependence for graphene is small below the Griineisen temperature Tbg given 
by Tbg ~ 5A^^/n^ K, where the concentration Ug is in unit of lO'^^cm"^ Jlii£. For our samples, 
Tbg ~ 54 K, and almost all our measurements are done below the Griineisen temperature. 

As kBTTtr/h ^ 1, we are by definition in the diffusive regime, for which EEI theory 
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FIG. 4: (color online) (a) phase coherence time extracted from the fit of the WL, see Fig. [3l 
Black circles: sample SI, red squares: sample S2. is roughly proportional to 1/T between 60K 
and 2.5 K and saturates at lower temperatures. The dashed line is a guide for eyes, (b) Coherence 
length versus 1/T, for the two samples. 

predicts a temperature dependent correction to the conductivity given by: 

AcTee = K-ee—r m (7) 

irn n 

where Kee is a prefactor whose value depends on the different channels contributing to the 
EEI^. In fig. [5|, we show 6axx(T) for samples 52 and SI. The expected log(T) dependence 
is evidenced and the slope gives Kee= 0.81±0.05 and 0.74 ±0.07 for samples 5*2 and 5*1 
respectively. These values are very similar to recent experimental findings on exfoliated 
graphene^i^. In a 2D system with a single valley, Kee takes the form:^ 

Kee = l + c(l - ln(l + F^)/F^) (8) 

where the unity represents the so-called 'charge' contribution (the Fock and singlet part of 
Hartree term), the c prefactor is the number of 'triplet' channels (from the Hartree term) 
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FIG. 5: (color online) experimental temperature dependence of C(T) = 6a'-' + 6a'^'^{T) vs log(T) 
for samples SI and 5*2. The slope gives Kgg = 0.81 and 0.74 for S2 and SI respectively. 

and Fq is the liquid Fermi constant. At low temperatures ksT < h/rgr (T < 15 K for 5*2 
and ?)QK for SI), when long range and short range scattering rates are important, the usual 
single- valley case is recovered with a prefactor c = 3. This leads to Fq = -0.13 ±0.04. This 
value is small and very close to the value of Fq recently found for exfoliated graphene^. 
This is somehow surprising, as the dielectric constant in SiC (10) is larger than in Si02 (3.9) 
and we would expect even smaller electron-electron interactions because of the screening 
of the substrate. However, numerical estimation of Fq following, for instance, Ref.- gives 
Fq ^ —0.09: a value compatible with our experiment. Also, Ref.— predicts values for 
the Fermi liquid constant only slightly smaller than our findings. Interestingly, it was also 
recently quoted^ that at intermediate temperatures h/rgr < ksT < h/rir, additional triplet 
channels originating from pseudospin conservation (ie valley degeneracy) become relevant 
and the number of triplet channels increases to c = 7. However, we do not observe any 
change in the slope in Fig. 5 from T=1.5K up to 50K. The situation is similar in Ref.—, 
where EEI is observed on a temperature range on which the slope should vary, but where 
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experimentally the slope remains constant. 

To conclude, we show, beyond classical geometric corrections, weak localization, weak 
antilocalization and electron-electron interactions in epitaxial graphene. Weak antilocaliza- 
tion is observable directly in the resistances, and its analysis gives access to the different 
scattering times, which are very close to those for exfoliated graphene on Si02 substrates. 
Electron-electron interaction gives also a small correction to the conductivity, which is not 
significantly smaller than in exfoliated graphene. 
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for a grant Juan de la Cierva. N. C. also acknowledges A. Bachtold's, A. Barreiro and J. 
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VI. APPENDIX A 

In Fig. [6], we present magnetoresistances of samples 5*1, 5*2 and of an additional monolayer 
53, more doped. A straight line aB has been subtracted for sample 5*2, for clarity. For all 
three samples, a dip is observed on the two sides of the central WL peak. We attribute the 
dip to WAL. 
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FIG. 6: transverse magnetoresistance for 3 different samples. A small depression attributed to 
WAL is visible around the WL peak for all three samples. 

VII. APPENDIX B 

We suggest that most part of the parasitic magnetoresistance comes from the invasive 
lateral probes. In order to sustain this assumption, classical magnetoresistances can be 
calculated by means of finite elements method (FEM) for the geometry of the devices. 
The simplest model only assumes different concentration and mobility under the lateral 
probes, and the results of this model are shown in Fig. [7] for sample 5*1. There is a good 
agreement with the experiment for a reduced mobility and an increased concentration under 
the probes, which seems reasonable if these region have been damaged during the electron 
beam lithography. 
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FIG. 7: (color online) experimental longitudinal and transverse magnetoresistances (only sym- 
metric or antisymmetric parts have been kept for clarity) for sample SI are compared to the 
result of a numerical model based on FEM. The model imposes an increased holes concentration 
= 7 • lO^'^cm"^) and a reduced mobility {fi = 0.06 m^V^^s^^) under the lateral probes. The 
rest of the graphene layer keeps its parameters as defined in table IB 
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